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ABSTRACT 

An X5.4 class flare occurred in active region (AR) NOAA11429 on 2012 March 
7. The flare was associated with very fast coronal mass ejection (CME) with its 
velocity of over 2500 km s” 1 . In the images taken with STEREO-B/COR1, 
a dome-like disturbance was seen to detach from expanding CME bubble and 
propagated further. A Type-II radio burst was also observed at the same time. 
On the other hand, in EUV images obtained by SDO/AIA, expanding dome¬ 
like structure and its foot print propagating to the north were observed. The 
foot print propagated with its average speed of about 670 km s -1 and hit a 
prominence located at the north pole and activated it. While the activation, 
the prominence was strongly brightened. On the basis of some observational 
evidence, we concluded that the foot print in AIA images and the ones in COR1 
images are the same, that is MHD fast mode shock front. With the help of a 
linear theory, the fast mode mach number of the coronal shock is estimated to 
be between 1.11 and 1.29 using the initial velocity of the activated prominence. 
Also, the plasma compression ratio of the shock is enhanced to be between 1.18 
and 2.11 in the prominence material, which we consider to be the reason of the 
strong brightening of the activated prominence. The applicability of linear theory 
to the shock problem is tested with nonlinear MHD simulation. 

Subject headings: magnetohydrodynamics(MHD) — Shock waves — Sun: corona 
— Sun: prominence — Sun: flares 


1. Introduction 

The solar corona is a place where plasma transient phenomena such as flares, coronal 
mass ejections (CMEs) and jets are continuously observed. These are vigorously investi- 
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gated in the framework of magnetohydrodynamics (e.g. ISh i bata fe Magaral (2 0111) 7 When 
magnetic energy is locally liberated in a short time scale as in flares or jets, magnetohy¬ 
drodynamic shock waves are generated. The shock waves sometimes propagate into the 
interplanetary space and reach the Earth. The arrival of the shock waves to the Earth some¬ 
times causes the sudden increase of radiation level and sudden change of Earth’s magnetic 
field by compressing it. 

The existence of coronal disturbances associated with flares are suggested by some obser¬ 
vation of far-away filament oscillation after flares. This is called ’winking fil ament’ after the 
obser vational feature where the filament appeared and disappeared in Ha ([Smith & Harvey 


19711 ). Sometimes, Type II solar radio bursts are observed associated with solar flares. Type 
II radio bursts are observed in dynamic radio spectra as bands slowly drifting down often 
in pairs differing in frequen cy by a factor o f abo ut 2. and interpreted as manifestation of 
coronal shock propagation ([Wil d et al. 19541: Ka i 1 9701 ). 


Globally propagating disturbances associated with solar flares are first observed in Ha. 
A fan-shaped ’wave fr onts’ were o bserved to propagate in the chromosphere after the oc¬ 
currence of a flare (IMoretom 119600 . This phenomenon was named ’Moreton wave’ after the 
name of the discoverer. The wave front of Moreton waves appear bright in line center and 
blue wing and dark in red wing followed by a fainter front whose signature of the intensity 
change is reversed (dark in the blue wing and bright in the red wing). This observational 
feature suggest down-up swing of chromospheric plasma during the passage of Moreton wave. 
Moreton wave propagated with speeds between 500 km s and 1500 km s -1 , and lasted for 
about ten minutes. In some events, Moreton wave front propagated faster than 2000 km s -1 . 
The sound velocity in the chromosphere is about 10 km s -1 . Fast mode magnetosonic wave 
speed is also the same order. If Moreton waves are magnetohydrodynamic wave propagating 
in the chromosphere, they have to be a strong shock wave with the Mach number of over 10. 
It was difficult to accept that such strong shock waves propagate a long distance without 
any significant dissipation. 


Uchida (119680 gave a comprehensive explanation of Moreton waves as a cross section 


between coronal MHD fast mode shock front and the chromosphere. This model of Moreton 
waves was widely accepted as ’Uchida’s Moreton wave model’ soon. According to Uchida 
model, there has to be a coronal counterpart of Moreton wave (namely coronal MHD fast 
mode shock front). People call the expected coronal counterpart of Moreton wave as ’invis¬ 
ible’ Moreton w ave. Winking filament were discussed as a result of the pass age of invisible 
Moreton waves ( Smith fc Harvev 1971 ; Eto et ah 2002 ; Okamoto et ah 2004 ). 


In 1990’s, wave-like disturbances propagating in the corona after flares were observed 
with the Extreme Ultra Violet (EUV) Imaging Telescope (EIT; Delabourdiniere et al. 1995) 
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on board Solar and Heliospheric Observatory (SOHO ; Domingo et al. 1995). This wave-like 
phenomenon was called 7 EIT wave’ after the name of the observational instrument. The EIT 
wave f r ont propagated concentrically from the flare site with the speed of about 250 km s _1 
( Thompson et al.l 1 9991 . 2000). At first, EIT waves were thought to be the expected coronal 
counterparts of Moreton wav es. However, EIT waves and Moreton waves showed qu ite dif¬ 
ferent observational features flEto et al.ll2002l; IShibata et al.ll2011[ IZhang et al.l 120111). which 


leads to a long standing discussion of the physical nature o f EIT wavesflWarmuth et al. 


2001; 

Vrsnak et a 

2002 

Gallaeher & Lone 

20 id) 


2002t iTripathi et al.l 120071 : IWarmuthl 120071: IWills-Davev fc Attrilll 120091 : 


Recently , with the Atmospheric Imaging Assembly (AIA; iTitle &; AIA team! (120061) : 
Lemen et al.l (120111) 1 on board the Solar Dynamic Observatory ( SDO\ Pesnell et al. 2012), 
fast coronal waves have been observed associated with flares (e.g. Liu et ah 20101 ; Ma et al. 


201l|). Chen & Wul (120111 ) found two different coexisting coronal waves, a slow coronal wave 
(i.e. EIT wa yel and a fast cor onal wave. These coronal disturbances are generally called 
EUV waves. Asai et al . (2 0121 ) reported a first simultaneous observation of Ha: Moreton 
wave, fast EUV wave and filament/prominence oscillations. In that paper, they showed that 
the fast EUV wave was the true coronal counterpart of the Ha Moreton wave. They also 
reported that after the Moreton wave disappeared, fast EUV wave propagated further and 
activated large amplitude filament/prominence oscillation. 

Large amplitude prominence oscillations caused by EUV waves are often observed by 
SDOj AIA. Such prominence oscillations are generally driven in a very short time scale and 
behave as damping oscillators, while in some cases, the oscillations lead to the eruption of 
the prominences. They are sometime s driven b y a su dden change of magnetic field structure 
through magnetic reconnections dVrsnak et al.l 120071 ). These oscillations have been investi¬ 
gated i n the context of (d i agnosis of physical q uantities and eruptive instability of the promi¬ 
nences (llsobe et al.ll2007l: ITripathi et al.ll2009f) . The prominence activations by coronal shock 
waves and its subsequent oscillation provide us with the insight into physical properties of 
both coronal shock waves and the prominences (such as the width of the wave front an d 
magnetic field structure supporting the prominences) (Gil bert et al. 120081: Liu et ah 2013). 


In this study, we analyzed the prominence activation process during the collision with 
an EUV wave associated with an X5.4 class flare on 2012 March 7. The unprecedented 
time and spatial resolution of multi-wavelength EUV data taken with SDO/AIA enables 
us to follow the prominence activation process in detail. We concluded, based on some 
observational evidence, that the EUV wave that activated the prominence is the same as the 
coronal shock wave propagating ahead of CME flank, which was seen in the images taken 
with the inner Coronagraph (COR1) of the Sun Earth Connection Corona and Heliospheric 
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Investigation (SECCHI) on board the Solar Terrestrial Relations Observatory ( STEREO ; 
Kaiser et al. 2008; Dreasman et al. 2008)-Behind satellite. We conld explain the compression 
and acceleration of activated prominence as a result of coronal shock wave transmission into 
the dense prominence. Furthermore, we estimated the coronal shock strength using initial 
velocity of activated prominence. 

In section 2, we overview the observation of X5.4 class flare and associated coronal 
disturbances. In section 3, we estimate the coronal shock strength using the initial velocity 
of activated prominence. 


2. Observation and analysis 
2.1. Overview of the flare and CME 

Figure 1 shows the soft X-ray light curve taken with GOES satellite. The first X5.4 
class flare occurred in active region (AR) NOAA11429 at 00:02UT and peaked at 00:24UT 
on 2012 March 7. One hour after the onset of the first flare, the second XI.3 class flare 
occurred at different part of the same AR. The two flares each were accompanied by fast 
coronal mass ejections (CMEs). The velocity of each CME is reported to be 2684 km s ” 1 
and 1825 km s -1 , respectively]^ The abrupt change of horizontal magnetic field and per¬ 
pendicular Lorentz force near the polarity inversion line (PIL) during X5.4 class flare are 
reported, and relationship betwe en the abr upt cha nge a nd the fast propagation velocity of 
the CME is pointed out f Petrie 2012 ; Wang et al. 2012h . In this study we investigate the 
coronal disturbances associated with the first flare of the two. 


2.2. A disturbance ahead of CME flank and Type II radio burst 

On March 7, the AR NOAA11429 was located at (N17,E16) in the heliographic coordi¬ 
nate. The STEREO-B satellite was then observing the Sun at 118 degrees east of the Earth 
having captured the CME just from the side. We analyzed the 5 minutes cadence data taken 
with STEREO-B/ SECCHI/COR1. 

Figure 2 shows the images taken with COR1 and EUVI on board the STEREO-B 
satellite. Figure 2(a) is the image taken with extreme ultraviolet imager (EUVI; Wiilser et 
al. 2004) 304 A pass band. In the EUVI 304 A image, we can see the polar prominence 


1 see, http://cdaw.gsfc.nasa.gov/CMEJist/UNIVERSAL/2012_03/univ2012_03.html 
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clearly. Figure 2(b) to (d) are the difference image of C0R1 image with embedded image 
of Figure 2(a). In Figure 2(b), a disturbance propagating ahead of expanding CME plasma 
appeared at 00:25UT. This disturbance reached the sky just over the polar prominence at 
00:31UT, and propagated further as shown in (c) and (d). The average speed of propagation 
of the disturbance from 00:25UT to 00:36UT was 420 km s _1 measured along the limb of 
the occulting disk of the COR1 images. 


On the other hand, in the dynamic spectrum in the range between 25MHz and 2500MHz 
taken with the Hiraiso Radio Spectrograph (HiRAS, Kondo et ah 1995), Type II radio burst 
was observed. The second harmonic component of the Type II radio burst appeared at 
00:19:12UT in 200MHz and disappeared at 00:30:36UT in 50MHz. From the average drift 


-l 


'Newkirk 

1961; 

Mann et al. 

1999) 


of the shock front ahead of the CME to be 672 km s 
Because of the coincidence, we regarded the coronal disturbance ahead of the CME bubble 
in the COR1 images to be the foot print of the shock wave which caused the Type II radio 
burst. 


2.3. EUV observation 

In order to investigate the time evolution of the coronal disturbances, we used the 
extreme ultra violet (EUV) images taken with SDO/ AIA. In this study, we used 2 pass band 
of AIA, namely 193 A and 304 A . The time cadence of AIA data is 12 seconds while we 
used 24 seconds cadence data in our analysis. 

Figure 3 shows the flare associated disturbances in 193 A images. At 00:18UT, dome¬ 
like structure appeared over the flaring AR. This dome-like structure kept on expanding and 
propagated upward from the solar disk. On the other hand, strong disturbance appeared on 
the north east of the AR at 00:23UT and propagated in the lower corona toward the north 
pole, and hit a polar prominence located there. By comparing COR1 and AIA data, we 
concluded that dome-like structure seen in AIA 193 A images is identical to the shock front 
ahead of CME seen in COR1 images. We consider the foot print propagated in the lower 
corona in AIA 193 A images are the foot print of the shock (Figure 4). 

These disturbances seen in EUV are both called ’EUV waves’. In the following, we 
distinguish the two disturbances as ’dome-like structure’ and ’foot print’ respectively. The 
dome-like structure expanded further and left away from the AIA field of view at 00:27UT. 
The foot print propagated further and hit the polar prominence at 00:33UT and caused 
large amplitude prominence oscillation. The average propagation velocity of the foot print 









6 


between 00:18UT and 00:33UT was 670km s 1 . We can also see that the wave front of the 
foot print more and more inclined to the propagation direction during its propagation. 

In Figure 3(b), the dome-like structure seems to already reach the polar prominence 
at 00:23UT. However, we notice that this is just because of a projection effect of three 
dimensionally expanding dome by comparing it with the shock front seen in STEREO- 
B/COR1 images. 


2.4. Prominence activation 

Taking a close look at the AIA 304 A images, we can clearly see how the large amplitude 
prominence oscillation was activated. Before the foot print ’hit’ the prominence at 00:33UT, 
the prominence was slowly swaying (Figure 5(a)). At 00:33UT, the foot print of the coronal 
shock wave hit the prominence (Figure 5(b)). When it hit the edge of the prominence, the 
prominence was strongly brightened. The brightness was about three times of that part 
before hit by the disturbance. After that the bright part propagated to the right in the 
image. The prominence was brightened from 00:33UT to 00:37UT, and during this time, the 
prominence material was accelerated from east to west in series (Figure 5(c)). 

We can also see that the prominence was accelerated perpendicular to the wave front in 
AIA 193 A images. In 193 A images, the prominence was seen as dark prominence (Figure 
6(a)-(d)). It was brightened just when the foot print arrived but not a long or strong 
brightening could be seen as in 304 A images. From AIA 193 A images, the propagation 
velocity of the wave front of the foot print just before the injection to the prominence was 
380km s _1 and the initial velocity of the activated prominence was 48 km s -1 . 


3. Estimation of shock strength in the corona 

We concluded that the foot print that activated the polar prominence was the shock 
wave which propagated ahead of the CME seen in STEREO-B/ COR1 images. Following are 
the reasons. 

First, the propagation location of foot print agreed well with the location of the shock 
wave ahead of CME bubble. The shock wave ahead of the CME appeared in STEREO- 
R/COR1 image at 00:26UT while the foot print in AIA 193 A images appeared at 00:23UT. 
The timing of its appearance agreed well within the time cadence of STEREO-B/ COR1. 
The shock wave in COR1 images passed the sky just over the polar prominence at 00:31UT 
while the foot print in AIA 193 A images hit the prominence at 00:33UT. We think that the 
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2 minutes delay is due to the inclination of the wave front towards the propagation direction 
in the corona. 


Second, the wave front of the foot print is inclined towards the propagation direction and 
the inclination increases with time.We thought this is because of wave refraction towards the 
photosphere because the phase velocity of MH D fast mode wave would increase with height 
in the lower corona ( Gopalswamv et al. 20011 1. The time evolution of the wave fron t agrees 
very well with the prediction of Uchida’s model of coronal fast mode shock wave flUchida 


19681 ). Because of the gravitational stratification, the plasma density of the lower corona is 


higher, which made the foot print of the dome-like coronal shock front bright in EUV pass 
band. 


Third, the polar prominence was strongly brightened, and at the same time, accelerated 
perpendicularly to the wave front. These behavior can be explained simply as a result of 
injection of compressive wave into the prominence. 

Looking at the prominence activation in the direction perpendicular to the wave front, 
we can regard this process as one dimensional fast mode shock wave transmission from corona 
into the dense prominence. In the linear case, the continuity of energy flux and momentum 
flux of the wave are written as 


Pc {V 2 - V 2 ) C c = Pp V t 2 C p (1) 

Pc(Vi-Vr)C c = Pp V t C p , (2) 


which reduces the equation 


v; 


1 

2 



V), 


(3) 


where p, V and C are plasma density, velocity amplitude of the wave, and phase velocity of 
MHD fast mode wave, respectively. The subscript i, t and r of V denote injected, transmitted 
and reflected waves, respectively, and the subscript c and p denote the quantities in the corona 
and prominence, respectively. 


Magnetic held strength around the prominence could be stronger than that of ambient 
corona, and therefore, the fast mode shock property in the corona could change gradually 
while approaching the prominence. If so, especially, wave refraction due to the spatial 
variation of fast mode propagation velocity in the corona around the prominence affects 
how strong the prominence is activated. In figure 7, we show the time-distance plot of the 
intensity of AIA 193 A images along the cutting-line shown in Figure 6 (a) from 00:00UT 
to 02:00UT. In Figure 7, on the other hand, we can see that the coronal fast mode shock 
front seen as bright signature approaching the dark filament propagates almost in constant 








speed. From this, we think that the spatial variation of magnetic held strength around the 
activated prominence did not affect much the process of prominence activation. Although 
the density discontinuously changes at the corona-prominence boundary, we think it is still 
reasonable to assume that the magnetic held strength at the corona-prominence boundary 
is continuous since the length scale of overall magnetic held structure is much larger than 
the width of the transition layer between corona and prominence. 

From the discussion above, we can express the compression ratio r c of the coronal shock 
wave using observed propagation velocity of injected wave C c and initial velocity of activated 
prominence Vt as 

r = = _1_ (4) 

c C c -v t l-ti + ys)^’ u 

where we express p p /p c as a. In order to simplify the problem, we treat the transmitted shock 
wave as a perpendicular one further on. If the local density gap between the prominence 
and corona is x and the volume filling factor of the prominence is fv, p P can be expressed 
using p c as p p = fvXPc + (1 ~ fv)pc- 

The relation between the fast mode Mach number Mf }C and the compression ratio r c in 
perpendicular shock case is 


M = i 2 r e ((2 — 7 ) r e + 7 (/3 C + 1 )) 

/,C ' ((7 + 1 )-^c(7- 1 ))(7/3 c + 2)' 


(5) 


In the linear problem, we also relate the compression ratio r p in the prominence with the 
one in the corona r c as 

2 y/cL 


r p = r c 


1 + 


( 6 ) 


We can see from equation ( 6 ) that the compression ratio of the shock wave is enhanced when 
it is transmitted into a dense prominence. 

From AIA 193 A images, we get C' c sin 6 = 380 km s _1 and Vjsin 0 = 48 km s ~ l where 
6 is the angle between wave propagation direction and the line of sight direction. Here we 
ass ume the local density gap x is 100 , the volume filling factor fv is between 0.001 and 
O.lf Terradas et al. 2008 : Labrosse et al. 2010 : Mackav et ah 2010 1. This leads to the range 
of the density gap p p /p c between 1.1 and 10.9. Estimated coronal shock properties are shown 
in Table 1. 


Especially, My c and r p falls into the range between 1.11 and 1.29, and between 1.18 
and 2.11 respectively. The brightness of the prominence material in AIA 304 A images is 
enhanced by a factor of about 2 during the activation, and it is roughly consistent with the 
result above if the brightening is due to compression. 
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4. 1.5 dimensional numerical simulation of fast mode wave transmission 

In order to estimate the coronal shock strength, we applied linear theory of wave trans¬ 
mission. In this section, we show the result of 1.5D MHD simulation of this model. 


4.1. Numerical settings 


We numerically solved 1.5 dimensional ideal MHD equations: 

dp _ _d_ ( T/ x 
dt dx {pVx) 


( 7 ) 


&v,) = -Iw + p + ±*1) 


( 8 ) 


3D y _ 3 (v o 

dt ~ dx (Vx y 


(9) 


3 A t/ 2 , P 
dt^2 pVx + 7 - 1 


+ ^ B ’ ) -~^i v A pv - + -,-i 


7 


1 

47T 


P + -B 2 ,) 


( 10 ) 


The numerical scheme is Harte n-Lax-van Leer-Discontinuities (HLLD) approximate 
Riemann solver flMivoshi fe Kusanol 20051) with second-order total variation diminishing 
(TVD) Monotonic Upstream-Centered Scheme for Conservation Laws (MUSCL) and sec¬ 
ond order Run ge-Kutta t ime integration with hyperbolic numerical divergence cleaning 
method flDedner et al.l 120021 ). The x-axis is uniformly discretized by N — 1000 grid points. 
Free boundary is applied as a boundary condition and the numerical box is x G [—1.0,1.0]. 
The ratio of specific heat 7 is assumed to be 5/3. The unit of length, velocity, time and 
density in the simulation are x c , Va c , t = x c /Va c and p c , respectively, where Va c is the Alfven 
speed. The unit of magnetic field strength is given as B c = a/ p c Uf c . The values of nor¬ 


malization units are set as a typical ones in the quiet sun corona: x c = 1.0 x 10 10 cm, 


Vac = 5.0 x 10 2 km s 1 , 

Be = 'JWac = 1-6 G. 


T = Xc/Vac = 2.0 x 10 2 


pc 


= 1.0 x 10” 15 


g cm 


-3 


and 


Initial conditions are as follows, 


P = 


p c = 1.0 (x < 0.5) 
p p (0.5 < x) 


( 11 ) 


P = P 0 


( 12 ) 
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B — (0, Bo) 

(2Vo sin(^), 0) (|x| < w) 

(0, 0) (|x| > w) 

where p c and p p denotes the density of the corona and prominence respectively 

Initially, gas pressure and magnetic field are set uniform. The magnetic field direction 
is perpendicular to the direction of wave propagation, since we think of perpendicular shock 
injection problem in this study. There is a density gap at x — 0.5, left and right side of 
which being the corona and the prominence, respectively. The density gap between the 
corona and the prominence is resolved by one numerical grid. The initial wave packet is the 
superposition of two sinusoidal-shaped wave packets propagating in the opposite directions 
to each other. We let the sinusoidal-shaped wave packet propagate in the corona and go into 
the dense region, ’the prominence’. 



(13) 

(14) 


For quasilinear cases Vo = 0.01, and for nonlinear cases Vo = 0.2. The width of the 
wave packet w is set to be 0.05. For each case, p c = 1, and p p ranges from 2 to 20. B c is 
fixed to be \/47r so that the initial Alfven velocity is unity in all cases. The initial pressure 
Po takes the value 0.025 or 0.1, each of which corresponds to plasma beta /3 = 0.05 and 
/3 = 0.2. The initial distribution of V x and p in quasilinear case (Vo = 0.01) with p p = 20 and 
is shown in Figure 8. While the wave injection into the prominence, some part of the wave 
energy is transmitted while some part are reflected (Figure 9). I 11 the nonlinear case, initially 
sinusoidal-shaped wave packet is deformed by nonlinear sharpening effect, and became shock 
before the injection into the prominence (Figure 10). 

In order to compare quasilinear and nonlinear cases, we define ’velocity transmittance’ 
of the wave T v as follows. 

T. = | (15) 

This is the ratio of the velocity amplitude of transmitted and injected waves (see Figure 
9 and 10). The velocity amplitudes of injected and transmitted waves are measured just 
before and after the injection, respectively. Linear analytic expression of this quantity as a 
function of the density gap a = p p /p c is (from equation (11)), 


T = 

- 1 - V 


2 

1 + y/a 


(16) 


Figure 10 compares the numerical results and analytic expression, (a) and (b) show the 
cases with (3 = 0.05 and /3 = 0.20, respectively. From the result shown in Figure 10, we can 
say that if the strength of the injected shock is not strong, the linear analytic treatment is 
not so bad. 
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Then, we define value R c as follows 

Re = (17) 

Po 

where po denotes the initial density distribution. At the wave front, R c corresponds to the 
wave compression ratio. The distribution of R c in quasilinear case just before and after the 
wave injection is shown in Figure 12, and the same plot for nonlinear case is shown in Figure 
13. The time evolution of velocity and compression ratio in quasilinear and nonlinear cases 
are shown in Figure 14 and 15. The steep overshoot at x — 0.5 after the wave injection 
is the result of the rightward shift of the prominence material, so it does not represent 
the wave compression ratio. From the linear theory, the compression ratio of the wave is 
expected to be enhanced after the transmission into the prominence (equation (6)). The 
expected enhancement of the compression ratio after the wave transmission can be seen 
both in quasilinear and nonlinear cases. 


5. Summary and conclusion 

In this paper we studied the nature of a globally propagated EUV wave associated with 
an X5.4 class flare occurred at AR11429 on March 7, 2012. The X5.4 class flare started at 
00:02UT and peaked at 00:24UT. A dome-like structure was observed in the images taken 
by SDO/AIA at 00:18UT and continued to expand. Then, at 00:23UT, another disturbance 
front appeared in the lower corona, and propagated toward the north with the velocity of 
around 670 km s -1 . The foot print reached the polar prominence at 00:33UT. On the other 
hand, in the images taken by STEREO/ COR1, a disturbance appeared at 00:25UT around 
the expanding CME bubble, and passed the north pole between 00:31UT and 00:36UT. The 
appearance of the disturbance agreed well with the start of detection of Type II radio burst 
by HIRAS. The time variation of the location of the EUV disturbance which propagated to 
the north also agreed well with the location of the disturbance seen in the STEREO/ COR1 
images. From these observational properties of the EUV waves, we consider the foot print 
as a footprints of shock front launched from expanding CME bubble. 

When the EUV wave hit the polar prominence, it is strongly brightened and the bright 
part propagated in the direction of EUV wave propagation. At the same time, the prominence 
was accelerated toward the direction of propagation of the wave. We consider the prominence 
activation as a result of shock injection. The strong brightening of the prominence could 
be explained as a result of the interaction of the shock with the prominence and resulting 
enhancement of the compression ratio of the prominence material. 

Using the observed velocity of activated prominence, the fast mode mach number Mf 
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of the coronal shock wave was estimated to be between 1.20 and 1.42 on the basis of linear 
theory. Also, linear theory predict the enhancement of the shock compression in the promi¬ 
nence material, which could explain the strong brightening of activated prominence in AIA 
304 A images. The estimated compression ratio of the shock in the prominence was between 
1.74 and 2.50. We checked the validity of the method by one dimensional numerical model 
calculation, and found that the linear theory is applicable when the shock is not strong. 


The SDO/AIA data are courtesy of NASA/SDO and AIA science team. The authors 
would also like to thank SMART team at Hida observatory of Kyoto university who con¬ 
tinuously provide high quality Ha full disk image of the sun. The simulation code used in 
this work is created with the help of HPCI Strategic Program. This work is financially sup¬ 
ported by JSPS KAKENHI Grant Numbers 24740331, 23340045 and 25287039. The autors 
are grateful to Dr. Andrew Hillier (Kyoto University) for helpful comments. 
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Fig. 1.— GOES Soft X ray light curve. X5.4 class flare occurred at 00:04UT and peaked at 
00:24UT. One hour later XI.3 flare occurred at different part of the same active region. 
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STEREO Behind EUVI 304 00:26:05 


STEREO Behind COR1 00:31:05 



STEREO Behind COR1 00:26:05 



STEREO Behind COR1 00:36:05 


(d) 


_ 71 

Dome 



Fig. 2.— Side view of the sun where we can see polar prominence and expanding dome-like 
structure. Images are taken by STEREO Behind telescope, (a): EUVI 304 A image. We 
can see the polar prominence clearly, (b)-(d): Difference image of COR1 image with EUVI 
304 A image embedded. Embedded 304 A image is the same as (a). The disturbance ahead 
of expanding CME appeared at 00:26UT and propagated to the north and passed the sky 
over the polar prominence. 
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Fig. 3.— (a)-(d): SDO/AIA 193 A difference image from 00:18UT to 00:36UT. (a): Dome¬ 
like structure appeared at 00:18UT. (b):foot print appeared at 00:23UT. (c): The foot print 
hit the polar prominence at 00:33UT. (d): foot print passed the polar prominence at 00:36UT. 
(e)AIA 304 image at 00:10UT. We can identify the polar prominence clearly, (f): Ha image 
taken by SMART in Hida observatory. We can see the Ha ribbons ob the flare and polar 
prominence. The white rectangles in (e) and (f) correspond to the FOV of Figure 5. 
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Fig. 4.— (a): SDO/AIA 193 AA image taken at 00:22UT. (b): Schematic figure of EUV 
disturbance seen in SDO/AIA 193 A image (a). The line of sight magnetic field obtained with 
SDO/HMI before the flare and the extrapolated coronal magnetic field (potential magnetic 
field) is shown. The red cross expresses the location of the prominence. 
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Fig. 5.— AIA 304 A image of prominence activation, (a): Before the arrival of the foot print 
at 00:31UT. (b): Just when the foot print reached the polar prominence at 00:34UT. The 
edge of the prominence strongly brightened, (c): While the acceleration of the prominence 
by the foot print at 00:36UT. A bright front propagates in the prominence toward the 
right. (d):When the prominence reached the maximum displacement at 00:55UT after the 
prominence activation had completed. The prominence kept oscillating after the activation 
had been completed. 
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Fig. 6.— AIA 193 A image of prominence activation. The FOV and timing of images (a)- 
(d) correspond to those of Figure 5 (a)-(d). The core of the polar prominence appears dark 
in 193 A pass band, (a): Before the arrival of the foot print at 00:31UT. (b): Just when 
the foot print reached the polar prominence at 00:34UT. (c): While the acceleration of the 
prominence by the foot print at 00:36UT. (d): When the prominence reached the maximum 
displacement at 00:55UT after the prominence activation had completed. The prominence 
kept oscillating after the activation had been completed. 
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Fig. 7.— AIA 193 A time-distance plot of prominence activation along the line A-B in Figure 
7. The pixels whose pixel count is lower than the threshold value which is set to emphasize 
both the signatures of wave front and dark filament is out of color table and appears white 
in the figure. The bright feature approaching the dark filament is the coronal shock front. 
We can see that the propagation speed of the shock front does not change much in the 
corona around the prominence. Just when the shock front arrives at the prominence, the 
prominence is suddenly accelerated. After the activation, the prominence reach its maximum 
displacement and continues oscillation. 
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Table 1: Estimated shock properties for 7 = 5/3 
on the basis of the linear theory 


/v a 


r C 

1 c 

'p 


0.001. 

0.05 

1.15 

1.18 

1.11 

0.001 

0.20 

1.15 

1.18 

1.11 

0.01 

0.05 

1.18 

1.38 

1.14 

0.01 

0.20 

1.18 

1.38 

1.13 

0.1 

0.05 

1.37 

2.11 

1.29 

0.1 

0.20 

1.37 

2.11 

1.29 


“The filling factor of the prominence 
& Plasma beta in the corona 

Compression ratio of the shock wave in the corona 
Compression ratio of the shock wave in the prominence 
“Fast mode Mach number of the shock wave in the corona 
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t = 0 


t = 0 




Fig. 8.— Initial profile of the density(left) and velocity(right) for a quasilinear case with 
P P = 20. 


t = 0.4 t = 0.8 




Fig. 9.— Vx just before and after the wave transmission in the quasilinear case. Velocity 
amplitudes of injected (Vj) and transmitted ( V t ) waves are indicated. 
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t = 0.35 t = 0.6 




Fig. 10.— Vx just before and after the wave transmission in the nonlinear case. Velocity 
amplitudes of injected (Vj) and transmitted ( V t ) waves are indicated. 
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Fig. 11.— The result of one dimensional simulation of shock transmission. The velocity 
transmittance T v is shown as a function of density gap a. Solid lines denote linear analytic 
solution, triangles and diamonds show the numerical results of quasilinear and nonlinear 
(shock) transmission cases, respectively, (a): Plasma beta /? = 0.05 case. Fast mode Mach 
number of the injected shock wave Mf = 1.16. (b): Plasma beta /3 = 0.20 case. Fast mode 
Mach number of the injected shock wave Mf = 1.17. 
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t = 0.4 t = 0.7 
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Fig. 12.— R c just before and after the wave transmission in the quasilinear case. Enhance¬ 
ment of R c after the transmission can be seen. 
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t = 0.35 t = 0.6 




x x 

Fig. 13.— R c just before and after the wave transmission in the nonlinear case. Enhancement 
of R c after the transmission can be seen. 
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Fig. 14.— Time evolution of V x and logR c in quasilinear case is shown. 
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Fig. 15.— Time evolution of V x and logi? c in nonlinear case is shown. 




























































































































































